Qubits as spectrometers of dephasing noise 
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We present a procedure for direct characterization of the dephasing noise acting on a single 
qubit by making repeated measurements of the qubit coherence under suitably chosen sequences of 
controls. We show that this allows a numerical reconstruction of the short time noise correlation 
function and that it can be combined with a series of measurements under free evolution to allow a 
characterization of the noise correlation function over many orders of magnitude range in timescale. 
We also make an analysis of the robustness and reliability of the estimated correlation functions. 
Application to a simple model of two uncorrelated noise fluctuators using decoupling pulse sequences 
shows that the approach provides a useful route for experimental characterization of dephasing noise 
and its statistical properties in a variety of condensed phase and atomic systems. 



I. INTRODUCTION 

A key step in the design of a quantum information pro- 
cessing device is to gain a quantitative understanding of 
the decoherence-inducing noise processes present in the 
system under study. Knowledge of the statistical proper- 
ties of this noise both informs and constrains theoretical 
models of the system, aiding in the design process. In 
some experimental qubit realizations, the technological 
implementation may allow for a direct measurement of 
the noise. In the superconducting flux qubit, for exam- 
ple, this is as simple as measuring the magnetic field fluc- 
tuations with a SQUID. Statistical properties of the noise 
may be readily computed from such a measurement and 
this has been very effectively used to measure the spec- 
trum of the ubiquitous 1// noise. Frequently, however, 
the noise acting on a quantum system is inaccessible to 
such direct measurements and the only recourse is in- 
stead to describe the noise indirectly through its effects 
on measurable quantities, such as a qubit decoherence 
rate. 

It has been recognized for some time that appropriately 
selected indirect measurements may provide insight into 
specific features of the noise source [HE]. For example, 
Schoelkopf et al. showed in 2002 that a single qubit is 
a valuable resource for measuring those characteristics 
of the power spectrum of an external noise source giv- 
ing rise to bit flips, i.e., dissipative characteristics of the 
noise source PQ. In many implementations, however, de- 
phasing is the dominant decoherence mechanism [3] and 
the techniques described in pQ are not directly applica- 
ble. Subsequent work has recognized the possibility of 
using pulsed spectroscopies as general diagnostic tools of 
spectral noise, including both dissipative and dephasing 
properties [2 . In particular, it has been shown that 
dephasing rates measured under application of dynami- 
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cal decoupling pulse sequences can allow differentiation 
between weak and strong environmental fluctuators [7]. 
However methods to make detailed characterization of 
dephasing noise that go beyond establishing such general 
features are still lacking. 

In this work we demonstrate that a controllable single 
qubit may be used as a sensitive spectrometer of dephas- 
ing noise that allows direct access to the noise correla- 
tion function. We first present a general procedure for 
estimating the correlation function of such noise using 
sequences of coherence measurements on a single qubit 
undergoing free evolution, and then show that the result- 
ing information may be extended to a signficantly larger 
range of time scales by use of pulse sequences that extend 
the coherence of a qubit by several orders of magnitude 

0- 

A major feature of our approach is the direct nature of 
the procedure. No fitting to parametric forms of corre- 
lation functions or spectral densities is required. Finally, 
we emphasize that, unlike the large literature of discus- 
sions of noise sources that focus on the form of the spec- 
tral density, our discussion here is focused on extraction 
of the noise correlation function. The present approach 
may be adapted to estimate the noise power spectral den- 
sity, as we describe in Appendix [C] However, since our 
reconstruction of the noise characteristics is explicitly lo- 
cal in the time-domain, our discussion is more naturally 
suited to the correlation function picture. 



II. MODEL 

We consider a single qubit subjected to a classical 
source of dephasing noise, described by the Hamiltonian 

H (t) = \a(t) ■ a + + V (t))a z . (1) 

Here a = (cr x , a y , a z ) are Pauli matrices and a(t) is a 
control field. For later convenience, we have separated 
the second term into a constant offset field, 770, and a 
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zero- mean stochastic process, rj(t). Such a Hamiltonian 
could arise, for instance, for the spin degree of freedom 
of an electron in a fluctuating magnetic field, a common 
source of environmental noise for many atomic qubit sys- 
tems, e.g., dopant spins in silicon [3]. We assume that 
the qubit can be initialized in an arbitrary pure state and 
that it can be measured in any basis. We will addition- 
ally assume the noise process to be wide-sense stationary, 
allowing us to write the correlation function as a function 
of only a single time, 

C{t) = ((m + v(t))(Vo + r?(0))) = Til + C n (t). (2) 

Here we have explicitly separated the full correlation 
function into its DC offset contribution, t]q, and zero- 
mean stochastic contribution, C n (t) — (r)(t)r)(0)) . 

Experimental constraints will necessarily introduce 
timescales that limit what may be learned about the sta- 
tistical properties of the noise from measurements on the 
qubit. The proposed experiment will consist of sequences 
of control pulses and measurements, each of which im- 
pose a specific timescale constraint. The longest of these 
timescales is the total length of time, T, over which the 
entire experiment is performed. No matter what mea- 
surements are made during T, correlations long compared 
to this time cannot affect these measurements and are 
therefore not accessible. A second constraint is that the 
state of the qubit must be reinitialized between each mea- 
surement, introducing a delay A t between experiments, 
or conversely, a ceiling on the maximum repetition rate 
to one measurement in time A t . The shortest time scale 
is that of the individual measurements. Since through 
the application of control pulses, the coherence of the 
qubit may be extended to a maximum time, T2, beyond 
which the coherence has decayed to a point where it is no 
longer measurably different from zero, we take this coher- 
ence time as the longest time available for a measurement 
of the qubit. 

Taken together, these limitations provide a natural 
separation of the problem into short timescales within 
the coherence time, t < T2, and long timescales lying in 
the range A t < t < T. Noise correlations occurring on 
timescales intermediate between T2 and A t are inacces- 
sible to measurement if T 2 < A t , as illustrated in Fig. [T] 
This puts limitations on the conversion of the directly 
determined noise correlations to a noise spectral density 
(see Appendix) and clearly implies that it is advanta- 
geous to make the measurement repetition time as short 
as possible. 

In the following section, we construct the experimental 
procedures necessary to extract the correlation function 
on i) times long compared to the measurement separation 
time, A t , and ii) times short compared to the maximum 
coherence time, T 2 . While both methods yield direct 
estimates of the correlation function, the two approaches 
differ in the type and sequence of both control pulses and 
qubit measurements that are employed. 
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Figure 1: (a) Illustration of the various timescales involved in 
measurement of the noise correlation function. The schematic 
shows a series of individual measurements (tall, black rectan- 
gles) distributed along a horizontal time axis. For simplicity, 
we have shown each measurement as lasting a time approxi- 
mately equal to T2 , the coherence time of the qubit. Qubit ini- 
tialization steps cause a delay At between measurements and 
the measurements are repeated until a final time, T. These 
times define the ranges which are accessible to direct measure- 
ment of the noise correlations (see text) and Tx denotes one of 
the timescales on which correlations are inaccessible to mea- 
surement, (b) Schematic reconstruction of a noise correlation 
function. For times shorter than T2, a continuous estimate 
is made by inversion of dynamical decoupling sequences (Sec- 
tion IV I. For long times, a discrete estimate is made at times 



equal to integer multiples of the measurement time, A t from 
measurements under free qubit evolution (Section III I. The 
time axis is scaled logarithmically here for greater clarity. In 
some systems, the short- and long-time estimates will overlap 
(i.e., T 2 > A t .) 



III. LONG-TIME CORRELATIONS 

To estimate the correlation function on timescales 
longer than the inter- measurement time, A t , we con- 
struct a sequence of qubit measurements following differ- 
ent times of free evolution, 5 t , according to the following 
protocol: 

1. initialize qubit in the +1 eigenstate of a x , 

2. allow qubit to evolve for a short time 5t, 

3. measure qubit in the basis of eigenstates of a y . 

We shall refer to these as free evolution (FE) measure- 
ments. We will assume that the j th run of this experi- 
ment begins at time, tj — jA t , i.e., the runs are equally 
spaced in time. This is a convenient but not essential 
requirement. During step (2) above, the qubit acquires 
a relative phase 
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(Vo + V(s)) ds. 
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Upon measurement in the a y 
probabilities are 



basis, the measurement dence on the control fields: 



i*(j) = r(lTsin« 



2 C 1 T <t>i) 



(3) 



Here we have taken a small angle approximation, which 
is valid if the total evolution time is chosen sufficiently 
small. The probabilities in Eq. ([3]) are approximately 
linear in the accumulated phase and may be used to es- 
timate the noise correlation function as follows. 

Making a large number, N, of repetitions of the above 
procedure will yield TV measurement results. We put 
these results in a vector, r, with rj = ±1 the result of the 
jth measurement. As shown in Appendix |aJ this mea- 
surement result vector may be used to estimate the full 
noise correlation function at times tk = kA t 



E[C(kA t )} = 
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The expected error of this estimate will scale inversely 
with number of measurements performed for each inter- 
val, 
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So for a given k, this error can be reduced by increasing 
N, the total number of FE measurements. 

Note that if we were to measure in the a x basis in- 
stead of a y , the measurement probabilities would then 
depend quadratically on the phase, P x {j) = cos 2 <j)j and 
P~(j) = sin 2 "(f)j. Such probabilities are independent of 
the sign of the acquired phase and are therefore not 
useful for determining the long correlations. However, 
repeated measurements in a x are nevertheless able to 
give useful information about the zero-time correlation, 
C(0) = (f? 2 ). In particular, by making a series of a x 
measurements and recording the values in a vector, r, 
we may obtain an estimate of the variance of 77, i.e., 
( V 2 ) = (l-{a x ))/(2df). 



H T (t) = ~( Vo +v(t))Ul(t)a z U a (t). 

Here U a (t) — e 1 ^ a ( s )-< jds / 2 j s the unitary operator deriv- 
ing from just the control field. The simplest dynamical 
decoupling procedures typically limit the control fields to 
7r-pulscs polarized along a x [10] . Since cr x a z <r x = —a z , 
the Hamiltonian in the interaction picture remains pro- 
portional to <j z and becomes 



H T (t) = \y{t) {va + v(t))vz 



(5) 



where we have introduced the pulse function, y(t), de- 
fined by 



y(t) 



1 after even number of 7r-pulses, 
— 1 after odd number of 7r-pulscs. 



Though we do not consider it here, the effects of nonzero 
pulse widths may be included to first-order through the 
modification, 

1 after even number of 7r-pulses, 
y(t) = ^ during application of 7r-pulses, 
— 1 after odd number of 7r-pulses. 

The pulse function describes the fact that, from the 
perspective of the qubit, each 7r-pulse acts to change the 
sign of the noise. When averaged over all possible noise 
trajectories, r](t), evolution under this Hamiltonian re- 
sults in dephasing of the qubit, which may be quantified 
by the decay of the expectation value of the coherence, 
(j + = (a x + ia y )/2. Taking the average over all possible 
noise trajectories yields 
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exp ( i I (770 + rj(s)) V( s ) ds ))Tr (cr+po) 



= exp Tr((7 +Po ) 



(6) 



IV. SHORT-TIME CORRELATIONS 

Noise correlations on timescales between t — and 
t = T2 cannot be investigated by the method described 
in the previous section, as these timescales are generally 
shorter than the minimum time between measurements 
(A t ). For these short timescales we show instead that 
judiciously constructed sequences of control pulses allow 
us to make direct measurement of the overlap integral of 
the noise correlation function with a filter function that 
is defined in terms of the applied control field. 

We begin by transforming the Hamiltonian into 
an interaction picture which removes the explicit depen- 



where we have made use of a cumulant expansion in the 
last line [TT] and the terms X (*) wn l be defined be- 
low. Note the two different sources of averaging for the 
qubit coherence. On the right hand side we have indi- 
cated averages over random variables, i.e., the stochas- 
tic average over all consistent trajectories of the noise 
term, with a single expectation value, (•). The double 
expectation, ((■)), of the coherence operator on the left 
hand side represents both this stochastic average over the 
noise realizations and the quantum average over the ini- 
tial qubit states, denoted by the usual Tr operation on 
the right hand side. For our purposes here, we are in- 
terested in experiments in which the qubit is initialized 
into the +1 eigenstate of a Xl a pulse sequence is applied, 
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and the coherence is measured at a time r. In this case, 
Tr((T + /?o) = 1/2. For notational convenience, we will 
drop the explicit dependance on t from the cumulant ex- 
pansion, assuming that it is always implicitly evaluated 
at t = r, at the end of the pulse sequence when the co- 
herence has maximally refocused. The first term in the 
cumulant series, the O th -order cumulant, is 



Correlation Filter Functions 



x (o) = _j r]0 



y(ti)dti 



This term, which is purely imaginary, represents the co- 
herent precession of the qubit due to the offset field. The 
0-order cumulant vanishes if we select any of the many 
refocusing pulse sequences for which J Q y(s)ds = 0. The 
next term in the expansion vanishes, 



= -i 



(v(h))y(ti)dt 1 =Q > 



because the stochastic term, rj{t), has zero mean by con- 
struction. In fact, all odd-order cumulants are purely 
imaginary and will vanish provided the unconditioned 
probability of a given noise realization, P(rj(t)), is sym- 
metric so that negative and positive contributions to the 
integral of y{t\) cancel. All even-order cumulants, I = 2 
and greater, are purely real and therefore contribute to 
decay of the coherence. These are therefore the terms 
that are responsible for dephasing of the qubit. The dom- 
inant decoherence causing term in the expansion ^ is 
the / = 2 cumulant 

x {2) (t)= [ dt x [ dh(v{tiHh))y{ti)y(t2). (7) 



Making use of the assumption of wide-sense stationar- 
ity that we imposed earlier, Eq. allows us to relate 
this expression to the stochastic part of the noise corre- 
lation function, C v (t2 — ti) = (r](ti)r](t2)) ■ By changing 
variables from ti,t 2 to u = t% — t\, v = t% + 1\, Eq. Q 
may be expressed as a single-variate integral over the 
stochastic part of the correlation function, 
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du C 7) (m)J : "(w), 



(8) 



with T(u) given as a second single-variate integral over 
a quadratic function of the control pulse sequence: 
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J ' [ " > = I dv v[ ) y ( 



(9) 



Eq. ^ is known as the coherence integral and defines 
the correlation filter function (CFF). The CFF specifies 
the regions of the correlation function that contribute 
to dephasing under a particular pulse sequence and is 
only defined for u e [0, r]. Some examples of correlation 
filter functions resulting from different pulse sequences 
are shown in Fig. [2j 




Figure 2: (Color online) Example correlation filter functions, 
i) (black, solid) Free evolution, ii) (red, dashed) Hahn echo, 
iii) (blue, dotted) 4 pulse Unrig [TDj sequence, UDD(4) 



We note that discussion of spin coherence decay is usu- 
ally given in terms of the noise spectral density (see e.g., 
|12|). which is the Fourier transform of the correlation 
function |13j . As explained earlier, the current formu- 
lation in terms of the time-domain correlation function 
is preferred here since it allows for consistency with the 
method of elucidating long-time correlations presented in 
the previous section. 

By taking the absolute value of Eq. |6j we can remove 
the dependance on the 0-order cumulant (since this is 
purely imaginary), resulting in 



l«°-+(r)»| 



exp -x 



(0) 



X 



(2, 



exp 



(-*») 



(10) 

This equation relates the CFF which is determined by 
and calculable from the pulse sequence, to the qubit co- 
herence and to the noise correlation function, Since the 
qubit coherence is measurable while the noise correlation 
function is unknown, this suggests that a direct estimate 
of the latter may be obtained from the former by a suit- 
able discretization and inversion of Eq. (10 1. We now 
show how this may be done numerically. 

To reconstruct the correlation function on short- 
timescales, one must measure the coherence decay for 
a number of different pulse sequences. Each pulse se- 
quence, i, determines a filter function Ti(t) as well as 
a coherence integral, xl = Jo* C^^Fi^dt. We define 
|((<7+))j| as the norm of the coherence measured after the 
pulse sequence associated with filter function Ti (t) . Now 
because the norm of the coherence is an experimentally 
accessible quantity, Eq 
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then allows us to calculate 
at time r. The procedure is 



the coherence integral \ 
then repeated for this particular pulse sequence 



suffi- 
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ciently many times to gather good statistics. By subse- 
quently repeating the entire procedure for a large number 
of different pulse sequences, {i}, we will gather a set of 
measured coherence integrals and filter functions from 
which an estimate of the short-time correlation function, 
E[C v (t)], may be regressed by making use of the theory 
of underdetermined least-squares |14j . This results in the 
expression 



t = T, 



£[^(t)]=£xfH^-(')' 



(ii) 



where F,-j = J °° Fi(t)Fj(t)dt is the filter function overlap 
matrix which may be constructed analytically for dynam- 
ical decoupling sequences or numerically for more general 
pulse sequences. This set of equations is under deter- 
mined because we are trying to reconstruct a continu- 
ous function by measuring a finite set of real numbers. 
Consequently, there are in general an infinite number of 
possible correlation functions which are capable of repro- 
ducing the measured coherences and it is necessary to 
impose an optimality constraint. Our estimate uses the 
Moore- Penrose pseudoinverse, F + , which yields the solu- 
tion with minimal Euclidean norm |14j . A derivation of 
Eq. [TT] is given in Appendix [D] 



Choosing pulse sequences 

The particular choice of pulse sequences will drasti- 
cally affect the quality of the correlation function esti- 
mate, both by dictating the range of time over which 
the correlation function may be measured and by influ- 
encing the accuracy of the estimate within that range. 
For instance, limiting oneself to a series of free evolution 
experiments will only allow for an estimate of the correla- 
tion function at very short times. Decoupling sequences 
greatly extend the average coherence time, facilitating 
a concomitant extension of the region of the correlation 
function that one can estimate with this procedure. 

To quantify the sensitivity of a particular set of pulse 
sequences to the correlation function at a particular time, 
consider a perturbation of the noise correlation function, 
Crj(t) —> C v (t) + X T 6(t — t). To first order, this pertur- 
bation changes the correlation function estimate to 

E[C v (t)] -> E[C n {t)\ + X T Y f J r i(t)Ff j ^ j (s). 

Taking the variation of this with A, gives a measure of 
the effect of the perturbation on the estimated correlation 
function: 



Squaring this quantity and integrating over t provides us 
with a positive scalar measure of the sensitivity of our 
reconstruction to variation of the correlation function at 



/>00 

Qto = £/ (WF+^it)) 2 dt (12) 
id Jo 

This quality function depends only on the filter functions, 
Fi(t), and the overlap matrix, F. Examination of the 
quality function for various sets of filter functions has 
empirically shown that the estimated correlation function 
becomes unreliable at times for which 



Q(t) < max(Q(t))/5. 



(13) 



V. NUMERICAL SIMULATIONS 

To illustrate the efficacy of our approach, we apply 
the procedure to a single qubit dephasing under the ac- 
tion of two mutually uncorrelated random telegraph (RT) 
fluctuators, using Monte Carlo techniques to simulate a 
statistically consistent noise trajectory for each measure- 
ment. Each fluctuator is capable of existing in cither of 
two states, ±??i, and will randomly transition from one 
state to the other at a rate 7, . To capture both the short- 
and long-time correlations that may be characterized by 
this spectrometry, we choose a fast, low amplitude fluc- 
tuator with parameters 771 = 1,71 = 10 and a slow, high 
amplitude fluctuator with parameters 772 = 10, 72 = 0.01. 
For simplicity, we set the offset field to zero, 770 = 0. 
The resulting noise correlation function can be calculated 
from ((BTJ to be C(f) = Tyfe" 2 ^*! + ^e~ 2 ^, as shown 
in Appendix [5] We simulated a series of N = 10, 000 
measurements, of which the first Ni — 5000 were free 
evolution decays to measure long-time correlations, while 
in the last iV s = 5000 we used dynamical decoupling to 
investigate the short-time correlations. The j th measure- 
ment was assumed to begin at time tj — j A t . For each of 
the Ni free evolution measurements, we simulated noise 
trajectories of length tpE = 0.04. For each of the N s dy- 
namical decoupling measurements a noise trajectory of 
length t n = 1 was simulated. The initial state of each 
simulated trajectory was conditioned on the final state 
of the previous trajectory according to Eq. (Bl I. 



th 



For the short-time correlations, the result of the 1 
measurement was simulated by first evolving the qubit 
under the combined action of a decoupling sequence and 
the i th simulated noise trajectory, then by randomly se- 
lecting a measurement outcome based on probabilities 
calculated from the usual Born rule. The pulse sequences 
chosen to investigate the short-time coherences are given 
in Table |TJ Because of their demonstrated success [T5] in 
extending coherence, we chose the Uhrig decoupling se- 
quence |10j as the basis for this simulated experiment. 

The resulting numerically reconstructed correlation 
function is shown in Fig. [3j where it is compared with the 
exact, analytical correlation function. The inset shows a 
plot of the corresponding quality function Q(t). Using 
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Pulse Sequence, Time Range, Divisions, Repetitions 



FE(1) 


0.1-0.5 


10 


100 


UDD(2) 


0.1-0.5 


10 


100 


UDD(3) 


0.1-0.6 


10 


100 


UDD(4) 


0.1-0.7 


10 


100 


UDD(5) 


0.1-0.9 


10 


100 



Table I: Pulse sequences used to measure the short-time cor- 
relations of the noise discussed in Sec. [V] For each pulse se- 
quence and each time step, one hundred noise trajectories are 
simulated and their effects on the qubit coherence is mea- 
sured. FE = free evolution, i.e., no pulses. UDD = Uhrig 
dynamical decoupling sequence [1UJ . 

the heuristic Eq. [13J we are able to disregard the short- 
time reconstruction of the noise for times log(t) > —1.3. 
We see that the reconstructed time correlation function 
of the dephasing noise demonstrates remarkable overlap 
with the analytic correlation function at both short and 
long times, validating the direct reconstruction approach. 



of donor spins and interfaces. 

Lastly, we note that the direct reconstruction method 
described here is not restricted to use with a single qubit 
and may also be used for measurement of dephasing noise 
acting on ensembles of qubits. In this situation however, 
unless the measurements can be spatially resolved, the 
reconstructed correlation function will be necessarily re- 
stricted to noise that is spatially correlated across the 
sample, such as that deriving from fluctuations in the 
applied magnetic field of a magnetic resonance experi- 
ment. 
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VI. DISCUSSION 



By inverting the conventional use of control pulse se- 
quences, we have shown that a single qubit can be a 
valuable resource for direct measurement of dephasing 
noise. In particular, we have demonstrated that direct 
reconstruction of short and long time noise correlation 
functions may be made using a combination of pulse se- 
quences and free evolution measurements. In this work 
we have employed dynamical decoupling pulse sequences, 
but the method could in principle be extended to more 
general sequences. 

Dephasing noise, pervasive in many quantum systems, 
is still relatively poorly understood. Our direct recon- 
struction method is general and not dependent on any 
specific physical features of the probe qubit. It may 
therefore be applied to any system for which dephas- 
ing is the dominant source of noise. It is particularly 
well suited to the measurement of dephasing noise at in- 
terfaces, e.g., for trapped ions or for dopants in semi- 
conductors. One significant possible application of this 
approach is to the measurement of interface noise ex- 
perienced by donor qubits in silicon-based devices. As 
shown in [3 , donor qubits in silicon near an oxide in- 
terface demonstrate a marked increase in coherence time 
as the distance from the oxide is increased. Theoretical 
models of the noise process causing this decoherence sug- 
gest that the presence of fluctuating dangling bonds at 
the interface is responsible for decoherence [9]. However, 
these theoretical models require a dangling bond density 
which is inconsistent with the measured density [§1 [TB] . 
The ability to make direct measurements of the statisti- 
cal properties of this noise could aid greatly in developing 
understanding of its microscopic origin and in construc- 
tion of new theoretical models to describe the interplay 
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Figure 3: Combined reconstruction of the correlation function of two mutually uncorrelated RT fluctuators at both short- 
and long-times, obtained with use of the set of pulse sequences given in Table [I] The dashed line is the analytical correlation 
function, the solid line is the short-time reconstruction and the black dots represent the long-time correlations. Inset is 
the quality function for the short-time reconstruction. Dashed-dotted lines demarcate low quality regions. The short-time 
reconstruction is unreliable at times for which the quality function Q(t) < 0.2, corresponding to log(t) > —1.3. This unreliable 
portion is separated by the dot-dashed line in the main figure. 



Appendix A: Estimate of long-time correlation 
function from free evolution measurements 



The correlation function at long-times is sampled by a 
series of free evolution measurements, yielding a results 
vector, r. We define the correlator of this result vector 



N-k 



(Al) 



The expected value of this correlator may be calculated 
from the measurement probabilities given in Eq. (§ ac- 
cording to 



i 



z=l m,n-±l 
N-k 

y 

N -k ^ 



4(iV — k) 



(A3) 



The covariance of the acquired phases may be simplified 

as 

{(t>i<t>i+k) = dt 1 dt 2 ((t} + r](ti)) (770 + r){t 2 ))) 



u 

2„2 



ti + k 

<ti+S t rti + k+St 

SNq+ I dti / dt 2 (r](t 1 )T](t2)) 
Ju Jti +k 

s tVo + / C n (u)f k (u)du. 

JkA-Si 



(A4) 



In the last equality we have changed variables in the in- 
tegral from t\ and t 2 to v — t 2 + t\ and u = t 2 — t\ and 
then integrated over v. The filter function, f{u), that 
appears in integral is defined as 

\/2(w - kA t + S t ) »e [fcA t - S t , kA t ] 
F k (u) = { V2(l-u + kA t ) ue[kA t ,kA t +5 t ] . 
otherwise 

However, from small 5 t , we can assume that C v (t) is con- 
stant over the range t € [kAt — $t, kA t + St] and so comes 
out of the integral. We can now rewrite the correlation 
function as 

{<f>i<fn+k) = $tVo + 6?C v (kA t ). 
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Combining this with Eq. A"3 we see that 

(C fc ) = 5?r 1 2 +S 2 C 1 (kA t ). 

Because the best estimate of (Ck) is the sample correla- 
tion, Cfc, given in Eq. |A1| we are left with 



N-k 



This may be solved for the full correlation function, 
C(kA t ), as 



^ N-k 



C{kA t ) ~ -2— — ^ 7 v r;r <+fe . 



<5 t 2 (JV-fc) 

In the very-long-time limit we expect that the stochastic 
part becomes completely uncorrclatcd, 

lim CM) = Um (r)(t)r](Q)) = 0. 

The stochastic part of the correlation function may thus 
be recovered from the full correlation function by sub- 
tracting, i.e., 



CJt) = C(t) - lim C(f). 

i— >oo 



which, for the transition rate matrix given above, simpli- 
fies to, 

C(t)=r) 2 e- 2 ^. 

Because the covariance vanishes, the correlation function 
of multiple, uncorrelated RT fluctuators is simply the 
sum of the correlation function for each fluctuator indi- 
vidually: 

CiAt) = (m(t)vj(o)) = (m(t)) (%(o)) = o. 

An alternative calculation for the random telegraph 
correlation function may be made as follows. In an in- 
finitesimal time, S t , the transition probabilities in the i th 
fluctuator are approximately linear in the rate: 

-Pflip = j8 t 

Pno flip = 1 - lS t . 

Therefore, the probability of n transitions in a time in- 
terval, A t is given as 



7 AA Ar ~ r V7Af 



AT! 



AT J nl(N-n)\ 



-7A 



, (7A t ) r 



Appendix B: Calculation of the correlation function 
of a random telegraph fluctuator 



The random telegraph fluctuator is defined as a clas- 
sical stochastic process taking one of two values, ±77, 
with a flipping rate, 7. The probability of being in the 
state ±77 at time t is given by p± (t) . Defining the vector 
p(t) = (p + (t),p-(t)), the time evolution of the probabil- 
ities may be written as 



dp(t) ( -7 7 



dt 



7 -7 



p(t) = rp\t). 



This equation defines the transition rate matrix, T. From 
this, the time evolution of the probability vector may be 
solved as 



p(t) = e t P(0) 



(Bl) 



The correlation function for the noise may then be cal- 
culated to be 

c(t) = ]T ViP(vi,tfa, o)vjP(vj, 0) 

U rt - 



The last term in the first expression above is a combina- 
torial factor. From this we see that the probability of an 
even number of transitions (which would leave the state 
unchanged) is 



P e (A t ) = Yl Pn(&t) = Ul + > 



-2 7 A t 



while the probability of an odd number of transitions is 



P (A t ) = Y Pn(^t) = \(1 
n odd 

The correlation function is then 



1,3 
,,2 



-(2P e {t)-2P (t)) 



Appendix C: Frequency domain filter functions 

The discussion in the main text for the short-time cor- 
relations may be instead expressed in terms of the spec- 
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tral density. Recall the I — 2 cumulant is 



/(2) 



V(ti)y(ti)dh / r)(t 2 )y(t 2 )dt 2 
o Jo 

dh / dt 2 (r)(h - t 2 )v(0)) y(ti)y(h). 
Jo 



Using the Wiener-Khintchine theorem, 

duj 



C(t) 



2n 



e lU3t S{u), 



we can rewrite the correlation function in terms of the 
power spectrum. This gives 



,(2) 



du> 

du) 
2^ 



S(u) 



e WT y{r)dT 



S(u})F{uj;t). 



(CI) 



The last line above defines the filter function as the 
square of the Fourier transformed pulse function, 



F(oj;t) = 



e lUT y(r) dr 



N 



1 + (-l) N+1 e iut + 2 ^(-l)V' 



>Ajt 



The filter functions indicate the range of frequencies of 
the noise power spectrum which contribute to dephas- 
ing. As mentioned in Sec. [TV] the effects of nonzero pulse 
widths may be included to first-order through the modi- 
fication 

y(t) = (_1)J?«(*')«B7* 

1 after even number of 7r-pulses, 
during application of 7r-pulses, 
-1 after odd number of 7r-pulses. 

Including this modification to the pulse function changes 
the filter function expression to 



F(w;t) = 



l + (-l) 



N 

,;=i 



j luiAjt 



COS ((jjtir/2) 



Appendix D: Estimation of correlation function by 
CFFs 



Suppose we have chosen a large number of pulse se- 
quences and constructed their associated filter functions, 
J-i(t). As described in the main text we have experimen- 



tal access to the coherence integrals, 

/>oo 

X f= / C v (t)F t (t)dt. 
Jo 

Because the filter functions are known in terms of the ap- 
plied pulse sequences, we can use this integral to describe 
correlation function, C v (t). If our set of filter functions 
were orthonormal, it would be a trivial task to expand 
the correlation function as a weighted sum of correla- 
tion filter functions, much like a Fourier series expansion. 
However, this is not the case, so we instead construct a 
new set of orthonormal functions. This can be done via 
the Gram-Schmidt orthogonalization procedure to yield 
the set 



(Dl) 



Properly normalized, these functions, fi(t), are orthogo- 
nal under the inner product 

/>oo 

(fi,fi) = / f i (t)f j (t)dt = S ij . 



We can now expand the stochastic part of correlation 
function, C^t), as 

i 

= X)cy<^(i) J C 1 ,(t))/ i (t) 

i,j 

i,j k 

= X T -c T c- f(t). (D2) 



Here, x and F are the vectors of measurement outcomes 
and filter functions, respectively, and c is the matrix of 



expansion coefficients from (Dl). We can determine the 
matrix c T c by examining the orthogonalized filter func- 
tions 

(/»(*)) = £ { c im^m{t),C jn J : n {t)) 

?7i, n 

= %■ 

These last two lines may be cast as a matrix equation, 
1 = cFc T . Here Fy = (J 7 ;, Fj) is the filter overlap matrix 
and 1 is the identity matrix. We see then that c T c = F -1 . 
However, F is likely to be numerically ill-conditioned, so 
we replace F _1 with F + , the Moore-Penrose pseudoin- 



verse of F. We can now rewrite (D2) as 



(D3) 
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In reality one will only be able to perform a finite num- 
ber of experiments, so this expansion is only approximate 
(much as a finite Fourier-expansion is only an approxima- 
tion of the expanded function) . We point out that we are 
able to express our estimate power spectrum entirely in 



terms of the filter functions, their overlap integrals, and 
experimentally accessible parameters. It is not necessary 
to explicitl y co nstruct the orthogonalized filter functions 
defined in (Dl), and we avoid the numerical errors com- 



monly associated [T7] with the Gram-Schmidt procedure. 
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